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A large variety of microscopic or mesoscopic models lead to generic results that accommodate 
naturally within Boltzmann-Gibbs statistical mechanics (based on Si = —k J dup{u) In p{u)). Sim- 
ilarly, other classes of models point toward nonextensive statistical mechanics (based on Sq = 
k[l — J du[p{u)]'']/[q—l], where the value of the entropic index g G K depends on the specific model) . 
We show here a family of models, with multiplicative noise, which belongs to the nonextensive class. 
More specifically, we consider Langevin equations of the type ii = f{u) + g{u)^{t) + ri{t), where ^{t) 
and ri{t) are independent zero-mean Gaussian white noises with respective amplitudes M and A. 
This leads to the Fokker-Planck equation dtP{u, t) = ~du[f{u)P{u, t)]+Mdu{g{u)du[g{u)P{u,t)]}+ 
AduuP{u,t). Whenever the deterministic drift is proportional to the noise induced one, i.e., 

/(it) = —Tg(u)g'{u), the stationary solution is shown to be P(u, oo) cc |l — (1 — q)f3[g{u)]'^^ 
(with q = and (5 — ) • This distribution is precisely the one optimizing Sq with the con- 

{J du [g{'a)\'^\P{u)\''} / du [P{u)Y^ = constant. We also introduce and discuss 



T + M 

straint {[g{u)\^)q 
various characterizations of the width of the distributions 
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Ubiquitous systems can be naturally described within 
Boltzmann-Gibbs (BG) statistical mechanics (based on 
the entropy 5*1 = —k J dup{u)hip{u)). These systems 
have in common the fact that they like equally to live ev- 
erywhere they are allowed to, that is, they are ergodic in 
the entire phase space. However, there are other systems 
that, depending on the initial conditions, may prefer a 
particular subspace. The BG scenario may not be appro- 
priate any longer and an extension of the usual thermo- 
statistical description, taking into account the features 
of such subspace, would be required. If that subspace 
has a scale invariant geometry, a hierarchical or multi- 
fractal structure, then the model points toward nonex- 
tensive statistical mechanics (based on the entropic form 
Sq = k[l- Jdu [p{u)]i]/[q - 1] , g e 3?) [1] (see [2] for 
reviews). Among the models which belong to this cat- 
egory, one finds low-dimensional dissipative and conser- 
vative maps [3], fractional and nonlinear Fokker-Planck 
equations [4], Langevin dynamics with fluctuating tem- 
perature [5], growth of many-body scale- free networks 
[6] and long-range many-body classical Hamiltonians [7] . 
The corresponding value of the entropic index q depends 
on the specific model, or, more precisely, on the nonex- 
tensivity universality class of the model. It is our purpose 
here to show a large family of models with multiplicative 
noise which belongs to the nonextensive class. 

Microscopic dynamics, containing multiplicative noise, 
may be encountered in many dynamical processes, such 
as in stochastic resonance [8] , noise induced phase transi- 
tions [9], granular packings [10], and others [11,12]. Due 
to its significance, stochastic processes with multiplica- 
tive noise have been the subject of numerous studies in 
the last decades [13-15]. Here we will consider processes 
subject to both additive and multiplicative noises and de- 
scribed by the dimensionless stochastic differential equa- 
tion of the form 



= f{u) + 9(u)m + Vit), 



(1) 



where u{t) is a stochastic variable, /, g are arbitrary 
functions (g{0) = 0), and ^{t),ri{t) are uncorrelated and 
Gaussian-distributed zero-mean white noises, hence sat- 
isfying 

mat'))^2MS{t^t'), {r^{tMt'))^2A5it-t'), (2) 

where M, A > are the noise amplitudes and stand for 
"multiplicative" and "additive", respectively. Clearly, 
some degree of correlation between the two noises can 
be of physical relevance, however, this remains out of 
the present scope. The deterministic drift f{u) can be 
interpreted either as a damping force (whenever u is a 
velocity-like quantity) or as an external force (when mo- 
tion is ovcrdamped and u represents a position coordi- 
nate). Other interpretations are possible as well, depend- 
ing on the particular system treated. 

It is interesting to note that additive and multiplica- 
tive terms in Eq. (1) could be gathered in an effective 
multiplicative noise term [16], however we prefer to keep 
track of both sources independently. 

The stochastic differential equation (1) is not com- 
pletely defined and must be complemented by an addi- 
tional rule. This is due to the fact that each pulse of the 
stochastic noise produces a jump in u, then the question 
arises: which is the value of u to be used in g{u)l This is 
the well known Ito-Stratonovich controversy [15,17]. In 
the Ito definition, the value before the pulse must be used, 
whereas in the Stratonovich definition the values before 
and after the pulse contribute in a symmetric way. If 
noise were purely additive, then both definitions agree. 

The Fokker-Planck equation for the probability den- 
sity P{u,t), associated to Eq. (1), can be ob- 
tained from the Kramers-Moyal expansion dtP = 
^„>;^(— 9„)"[r''"^P], where the coefficients are given by 

t) - -^linv^o^^^^^^^^ These coeffi- 

cients can be readily obtained following the standard 
lines found for instance in [15]. Using the Stratonovich 
definition of stochastic integral, one gets 
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D(^^\u,t) ^ f{u) + Mg{u)g'{u) = J{u), (3) 
D(^\u,t) = A + M[g{u)]^ = D{u), (4) 

while (u, t) = for n > 3. Then, one arrives 
straightforwardly at 



dfP = -duj{u). 



(5) 



where j{u) = J{u)P — du[D{u)P] is the current. In the 
Ito calculation, Eq. (3) becomes J{u) = f{u), that is, the 
noise-induced (or spurious) drift is missing. In what fol- 
lows we will adopt the Stratonovich definition. However, 
this choice will not affect the present discussion except- 
ing for a redefinition of some of the parameters that are 
involved. 

Eq. (5) can also be written as 

dtP = -du{f{u)P) + Mdu{g{u)dulg{u)P]) + Ad^uP 

(6) 

In some processes, the deterministic and noise-induced 
drifts may have the same functional form. Let us put 
this condition as follows: 



/(u) = -Tg{u)g'{u), 



(7) 



T being a proportionality constant. In other words, f{u) 
is derived from a potential-like function V{u) = ^[glu)]"^. 
Let us note that the particular case g{u) oc f{u) cx u, 
which is a natural first choice for a physical system, veri- 
fies this condition. However, since no extra calculational 
difficulties emerge, we will discuss here the more general 
case (7). Notice that in the absence of deterministic forc- 
ing, condition (7) is trivially satisfied for any g by settig 
T =0. 

In the present paper we will restrict to the station- 
ary solutions for no flux boimdary conditions (i.e., such 
that j{—oo) = j{oo) = j{u) = 0), although more general 
conditions could in principle also be considered. If (7) 
is verified, then, the stationary solution Ps{u) is of the 
^'-exponential form appearing in nonextensive statistical 
mechanics [1]. More precisely, in that case, one obtains 



Ps{u) cc l + {q-l)l3[g{uf 



(8) 



where (3 = = '^^^ (where by T we generically mean 
the amplitude of an effective noise) and 



q = 



T + 3M 

T + M 



(9) 



As an asid(^ comment let us mention that if we had used 

the Ito convention we would have obtained q = ^^-r44f ■ 

Let us now go back to our Stratonovich choice. 

Ps is normalizable if \g{u)\ grows with |u| faster than 
1 

|y|i+T/M (^7- > —M). This probability distribution func- 
tion (pdf) optimizes 



with the constraint 

_jdu[g{u)r[p{u)r_ 



Jdu[P{u)]^ 



(10) 



constant. (11) 



The condition (7) is not necessary for having solutions 
of the (/—exponential form, as it follows along the lines 
of [18,19], in spite of the fact that the models therein 
considered are different from the present one. 

Let us consider in more detail the case when both 
the forcing and the multiplicative noise depend on the 
stochastic variable u as a power law, that is, when the 
Langevin Eq. (1) becomes 



ii = -7«|u|'-^ + u\u\^-'^^{t) + r]{t), 



(12) 



with r, s > and a drift coefficient 7 typically posi- 
tive. In this case, the deterministic drift is derived from 
a confining potential-like function of the form V{u) = 
7|u|''+^/(r -|- 1). The corresponding Fokker-Planck equa- 
tion becomes 

dtP = -du{-7u\ur^P) + Mdu{\u\'du[\urP]) + AduuP 



Hence, its stationary solution Ps{u) is given by 

Po e~'*("^ 

Po being the normalization constant and 



(13) 



(14) 



h{u) = 



7|u 



l+r 



A[l + r] 



2s 



r l+r 



2s 



(15) 



where 2F1 is the hypergeometric function. 
We shall analyze now some limiting cases: 



(A) For vanishing deterministic forcing (7 — » 0), Eq. (14) 
becomes 



Ps{u) = 



with 



Pn = 



(l+M|^|2.)l/2 

(M/>i)i/(2-) r(i) 
2r(i + i)r(i-^) ' 



for s > 1. 



(16) 



(17) 



This pdf is of the ^-exponential form (with g = 3) as 
expected because, in this case, condition (7) is trivially 
true since it corresponds to r = 0. Interestingly, in the 
presence of multiplicative noise, the steady state Pg is 
normalizable even in the absence of a confining potential, 
as long as s > 1. The so called spurious drift which origi- 
nates from the multiplicative noise term in the Langevin 
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equation is responsible for kicking the system back close 
to the origin, where fluctuations are smaller. On the 
other hand, the additive noise plays an essential role, 
providing fluctuations which avoid the full concentration 
(at the origin) that would occur otherwise. This type of 
stabilizing effect of the multiplicative noise is long known 
[20]. 

(B) In the limit M — !■ 0, all other parameters being fixed, 
i.e., when no multipliciative noise is present, Eq. (14) be- 
comes of the following stretched exponential form: 



Pn = 



T(- + — !— ) 



2Ar(i + J-)r( 



7 
2sM 



(22) 



Psiu) = 



2r(i 



In particular, for the linear forcing (r 
pdf is recovered. 



(C) In the limit ^ — » 0, i.e., for vanishing strength of the 
additive noise, the steady state is normalizable for s < 1 
and p = r+ l — 2s>0. The condition p > implies 
that the potential of the drift has to be steep enough to 
confine the system and yield a steady state. Then, we 
have 



Ps{ii.) is normalizable for s + -y/M > 1. Notice that q- 
exponcutial pdf's can appear even for negative 7 (i.e., 
for repulsive deterministic forces). The 1 + r = 2s class 
includes the particular linear case r = s = 1, that has 
already been treated in the literature [11,21]. In Fig. 1 
we exhibit the steady state pdf (20) for several values of 
the system parameters. 

Assuming that (u) = 0, the width of P{u,t) can be 
characterized in many ways, such as (i) the inverse of the 
height at the origin, namely 1/P{0,t); (ii) the width at 
half height. A; (iii) the square root of the mean value of 
u^; (iv) the 2s-root of the mean value of (v) the 

1) the Gaussian square root of the g-expectation values of u^; (vi) the 
2s-root of the g-cxpcctation value of or even com- 

binations of these, for instance (vii) ■y/A/P(0, t) (see also 
[22]). Let us recall that the normalized g-expectation of 
ip is defined as 



(18) 



^^^_ fdui^iu)[P{u,t)]'i 



(23) 



Psiu) 



P [life]' 



2r(i^) 



(19) 



Vanishing A concentrates the probability at the origin. 
Again, the Gaussian distribution is recovered for r = 1 
(harmonic forcing) and s = (additive noise). 

The limits analyzed above are not generically inter- 
changeable, that is, convergence is not necessarily uni- 
form. 

In the general case where {A, M, 7) are all finite, the 
same dependcmcc on u as that in Eq. (19) is obtained 
for sufficiently large \u\. In fact, if one defines the dimen- 
sionless variable u = u/X with A = {A/M)^^ ^ it is clear 
that the asymptotic expression of Eq. (14) for \u\ — > 00 
corresponds to both limits juj — > 00 and ^ ^ 0. 

For finite {A, M, 7), in the p = marginal case where 
the drift and multiplicative-noise exponents are related 
through 1 -h r = 2s (hence condition (7) is verified with 
r = 7/s), a g-exponential pdf emerges at asymptotically 
long time. In fact, in that case, the hypergeometric in 
Eq. (15) becomes 2-Fi(l, 1; 2; —z) = ln(l -|- z)/z and one 
gets 



hence, {1(^)1 equals the usual mean value (ip). 

For the usTial case where P{u, t) is a Gaussian, all these 
definitions basically coincide. This is not so in general, 
as wc shall illustrate in what follows for the stationary 
state of Eq. (20). 

(i) The inverse of the height at the origin: 



A 



2r(i + J-)r(i + ^-J-) 



(ii) The width at half height: 

A = 2A {2^/^27M+i) - 1)17. 

(iii) The square root of the mean value of u^: 



i^/^ = A 



^(is) + 2sM ^ 2s) 



(24) 



(25) 



(26) 



for 7/M > 3 — s (otherwise it diverges), 
(iv) The 2s-root of the mean value of [mI^*, 



{\nn 



2s\l/2s _ 



A 



(7/M - s - 1)- 



(27) 



Psiu) 



Po 



(1 + |«/A|2-)^ 



with 



and 



q = 



(7/s) + 3M 
(7/s) + M 



(20) 



(21) 



for 7/M > 1 -|- s (otherwise it diverges) . 

(v) The square root of the ^-expectation value of u^: 



\2s) ^ \2 ^ 2sM 2s > 
^\2s' V 2 ' 2sM 2s > 



(28) 



for 7/M > 3(1 — s) (otherwise it diverges), and 
(vi) The 2s-root of the g-expectation value of \u\'^^ 



3 



{\uW = x> (29) 

(7/M + S-IY' 

for 7/M > 1 — ,s. Notice that this condition has al- 
ready been encountered. Indeed, it is necessary for nor- 
malizability. An important remark is mandatory: of all 
the above characterizations of width which arc based on 
expectation values (i.e., (w^)^/^, (Iwp**)^/^**, (u^)^^ and 

), only the last one does not diverge in any 
admissible (i.e., normMizable) case. This is particularly 
remarkable because this generalized expectation value is, 
as already seen, precisely the one to be used as constraint 
in the optimization of Sq. In addition to this, it is wor- 
thy noticing that (Iwp'')^^* is in all cases comparable to 
■\/A/Ps(0), a semi-empirical quantity which takes into 
account the contributions of both body and tails. 

In Fig. 2 we exhibit all these quantities as a function 
of M/7 for typical values of s with r = 2s — 1. As M/7 
increases, in all cases, Ps{u) becomes more peaky around 
the origin (A decreases). Moreover, for increasing M/7 
and s = 1, the tail tends to l/\u\ (hence Ps(0) — > 0). 
Also for s = 1/2, Ps(0) ^ as the pdf becomes more 
tailed, but M/7 has the upper bound 2. Contrarily to 
the previous cases, for s > 1, as M/7 increases, Ps(0) 
also increases. 

The g-exponential character is not exclusive of the 
steady state but it also emerges along the time evolu- 
tion of the pdf, as illustrated in Fig. 3. In this figure 
we employ the semi-ln^ representation, where In^x = 
{x'^-t - l)/{q - 1). Notice in Figs. 3.b and 3.c that the 
curve for f ^ 00 is a straight line. In 3.c, curves are 
almost straight lines, thus indicating that the pdf's are 
very close to (/-exponentials along time. 

In the presence of multiplicative noise, the system vari- 
ables directly couple to noise. Therefore, behaviors are 
observed that can not occur in the presence of additive 
noise alone. In particular, g-exponential pdf's of the form 
given by Eq. (8) can arise. For the present processes, q- 
exponentials do occur either in the absence of forcing 
(7 = 0) or for drifts verifying Eq. (7), which, for the 
power-law case, becomes r = 2s — 1. The determinis- 
tic forcing does not need to be confining (i.e., 7 < is 
allowed) for the formation of a (/-exponential stationary 
pdf, provided 7/M > 1 — s. In any case, if a repulsive 
effective force prevails for some period of time, the vari- 
able u can take large values and power law tails can arise. 
In particular, a ^'-exponential occurs when the stochastic 
forcing is proportional to the deterministic one. Alter- 
natively, during the intervals when the effective force is 
attractive, the probability tends to concentrate at the 
origin. Then, at this stage, the additive noise plays a 
fundamental role allowing the existence of a normalizable 
steady state by avoiding collapse of the pdf at the origin. 
When both kinds of noise occur simultaneously, the pres- 
ence of multiplicative noise can not be formally avoided 
by a simple transformation of variables. The particular 
interplay between additive and multiplicative noises as 



well as that between deterministic and stochastic drifts 
can lead to the appearance of g-exponentials. 

It is worthy to emphasize at this point that the sta- 
tionary solutions of the present problem have the general 
form given by Eq. (14). The (/-exponential pdf's repre- 
sent a special case, which in turn includes the Boltzmann- 
Gibss pdf as an even more special one, corresponding to 
the standard thermal equilibrium. Furthermore, the q- 
exponential pdf's are unique in the sense that they op- 
timize, under appropriate constraints, the entropy func- 
tional Sq (Eq. (10)). This entropic form is in turn unique 
in the sense that it is the only one which satisfies [23] a 
set of conditions naturally generalizing both the Shan- 
non and the Kinchin axioms. Also, Sq is consistent with 
stability (or robustness) of the (/-exponential pdf's in the 
sense described by Abe [24], whereas Renyi entropy and 
the normalized version of Sg are not. 

Let us conclude by stressing that (/-exponentials have 
also been observed in a variety of similar processes 
[5,18,25]. The mechanism leading to such distributions is 
expected to be present in systems with long-range mem- 
ory, long-range interactions, fractal or hierarchical struc- 
tures and similar scenarios. 

In what concerns the model defined by Eqs. (3)- (5), 
it is clear that f(u) and g{u) are generically indepen- 
dent functions. However, it does occur that when they 
are connected through Eq. (7), the g-exponential form 
emerges naturally. It would no doubt be very interesting 
if we had a geometrical interpretation of this fact. Hints 
along this line would be welcome. 

As illustrated in Fig. 3 (deterministic drift propor- 
tional to the stochastic one), tlic^ (/-exponential is exact 
for t —^ 00 and it is an excellent approximation \/t. 

The best characterizations of the width of the pdf's 
clearly are those which remain finite under generic cir- 
cumstances. In our case, this happens for 1/P(0,i), A, 
(1^1^*)^^** or combinations such as y^A/P(0, t). The lat- 
ter two are preferable since they contain attributes of 
both body and tails of the pdf. For a mathematically con- 
venient and generic characterization, clearly (|Mp'')J^'^* is 
the most appropriate. 
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FIGURE CAPTIONS 

Fig. 1. Steady state pdfs for s=0.5 (a), 1 (b), and 2 (c), 
with r = 2s — 1, A/^ = 1 and different values of M/7 
indicated in the figure. Notice that the alternative rep- 
resentation XPs{u) vs. u/X with A = (A/M)^/^" would 
give profiles independent on the particular choice of A/'y. 

Fig. 2. Different parameters characterizing the "width" 
of the distribution as a function of M/7 for s=0.5 (a), 1 
(b), and 2 (c), with r = 2.s - 1 and A/7 = 1. For 7 > 0, 
normalizabity implies M/j < 1/(1 — s) if s < 1; for s > 1, 
the physical region of M/7 extends up to infinity. 

Fig. 3. Time evolution of the pdf for r = s = 1, A = 

7=1 and M = 0.2 at different times indicated in the fig- 
ure, in different representations: linear (a), semi-ln^ with 
Q = Ty^^M (^)' ^^'^ semi-ln,', with time dependent q' (c). 
We used the ansatz P = Po{l + {q' - l)/Ju2)i/(i-g'), ^ith 

^° ^ V /•j(g'^-i) p^^ 'i ^ - ■ In the inset we present q' and 

13 vs. t. In particular, ^(0) = 1, l//?(0) = 0, q{oo) = 4/3 
and l//3(oo) = 5/3. 
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